Astronomy Sz Astrophysics manuscript no. 0458-astroph 


January 20, 2009 


(DOI: will be inserted by hand later) 





Why did Comet 17P/Holmes burst out? 
Nucleus splitting or delayed sublimation? 
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p-^ ' Abstract. Based on millimeter-wavelength continuum observations we suggest that the recent "spectacle" of 

comet 17P/Holmes can be explained by a thick, air-tight dust cover and the effects of H2O sublimation, which 

i-C ■ started when the comet arrived at the heliocentric distance < 2.5 AU. The porous structure inside the nucleus 

provided enough surface for additional sublimation, which eventually led to the break up of the dust cover 
and to the observed outburst. The magnitude of the particle burst can be explained by the energy provided 
by insolation, stored in the dust cover and the nucleus within the months before the outburst: the subliming 
surface within the nucleus is more than one order of magnitude larger than the geometric surface of the nucleus 
- possibly an indication of the latter's porous structure. Another surprise is that the abundance ratios of several 
molecular species with respect to H2O are variable. During this apparition, comet Holmes lost about 3% of its 
I , mass, corresponding to a "dirty ice" layer of 20 m. 

^\ ' Key words, comets: general - comets: individual: 17P/Holmes 

m \ 
t> ■ 

. 1. Introduction Montalto et al. (|2008p report a significant disassembly 

of the nucleus, not even excluding a complete disintegra- 

Comet 17P/Holmes was serendipitously discovered during t ion. Earlier, Sekanina ((15521 had classified types of split- 

an outbreak on 18 92 November 6 by Holmes (reported ting comets: (a) single comets that break up into two or 

by Plummer [1893J while he was observing the nearby mor6j ( b ) come ts that disintegrate or suddenly disappear, 

Andromeda galaxy (M31). Until early 1893 January, the ( c ) an d those with a pancake-shaped companion nucleus 

comet faded from magnitude 4 to 9-10, after which a that disintegrates into microscopic dust grains. Recently, 

jJJ ■ second eruption to w 5 mag occurred. Obviously, this Sekanina {Ugl summarized the optical observations of 

C3 ' light curve is different in time dependence and amplitude 17P /Holmes and some other comets for comparison. All 

from what was observed during the most recent appari- types of splitting comets start in his hypothesis with a 

tion (2007/8). After the early observations, Comet Holmes ma j or outburst. The "megaburst" of 17P/Holmes is of 

was lost for some years, but later on recovered as an al- type (c), starting with an exothermic reaction, resulting 

most "dead" comet with magn itudes of w 16-17 near ev- m a rap idly expanding cloud of microscopic dust parti- 

ery perihelion. Whipple ((19860 analyzed the historic data c i es . But not all major outbursts end in splitting: e.g., the 

again and explained the two outbursts by grazing encoun- one of comet Halley on 1991 February 12 at a heliocentric 

ters of a small hypothetical satellite with the nucleus: the distance, r, of 14 AU (Sekanina et al. [T9921 . 

first one on 1892 Nov. 4.6, 1892 and the second on 1893 

January 16.3.; even though these encounters could not be ~ ■ 
„ , . . r , , . . , 2. The nucleus 

confirmed, his review of the historic observations allows 

this event to be discussed again in connection with the 2.1. Time line 
latest outburst discussed in this paper. 



o 



The "engine" behind the cometary activity of Comet 

Holmes is the production of gaseous water as described 

W. J. Altenhoff by Delsemme (1982). Its production rate, Q(H 2 0), is a 
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function of heliocentric distance. It is oc 1/r 2 for low val- 
ues of r, while for r > 1.5 AU the dependence becomes 
highly nonlinear. Dclsemmc defines the limit of sublima- 
tion r , the heliocentric distance beyond which 97.5% of 
the energy received by insolation is re-radiated, and only 
< 2.5 % are used for vaporization. For water ice, r$ is 
about 2.5 AU. 

Comet Holmes is a short-period comet in the Jupiter 
family (a JFC). Its average perihelion distance, q, over the 
last 6 apparitions was ss 2.2 AU (Marsden and Williams 
1999}, close to the limit of H 2 sublimation. But the 
perihelion distance of the most recent apparition was 
at a q of 2.05 AU. With the steepened production rate, 
mentioned above, the H2O production is increased by a 
factor of 2. This is possibly responsible for the outburst. 

The outburst happened according to Hsieh et al. 
(j2007j) on 2007 October 23.8, 173 days after perihelion pas- 
sage, or 361 days after crossing r = 2.5 AU. Probably - be- 
cause of the low-level cometary activity - the nuclear sur- 
face was free of ice and the icy nucleus was covered by some 
sort of a rubble pile (Jewitt, [1992) or dust-particle man- 
tle, causing the delay of visible cometary sublimation by 
months. During this period, the dust cover was "air tight" , 
preventing the sublimated gas to escape. Sublimation in- 
side the nucleus continued until the gas set free by this 
process broke up the dust mantle - the "outburst" . 

2.2. Model parameter 

Dust cover. All accurately measured cometary nuclei show 
a geometric albedo, p between 0.02 and 0.05 (see e.g. 
Jewitt l2005|) . which is indicative of a dust cover. If closely 
packed, this material is thought to have a density, p, rj 
1 g cm" 3 , as inferred from numerous radar observations 
(Harmon[l999; Harmon et al. 2005). Information from the 
collision caused by the Deep Impact mission revealed that 
the nucleus of comet Tempel 1 had a devolatilized dust 
cover of about 1 m, with very little H2O inside and none 
on the outside (Sunshine et al. I2007P - identical to what 
we assume for 17P/Holmes. Below the dust cover of comet 
Tempel 1, an at least 10 m thick layer of fine grained wa- 
ter ice particles was found, which appeared to be free of 
refractory impurities! This "clean" ice may originate from 
repeated sublimation and deposition inside the enclosed 
nucleus approaching and leaving the solar neighborhood. 
It is likely that comet Holmes has a similar layer. But for 
our model we neglect this detail and assume for the mass 
estimates "dirty" amorphous H2O ice throughout the nu- 
cleus. 

Diameter. Until recently, the resolution of optical tele- 
scopes was not high enough to directly measure the nu- 
clear diameters of JFCs. Instead, absolute magnitudes of 
the nuclei were determined and nuclear diameters were 
calculated, assuming a geometric albedo p = 0.04, because 
observations constrain the albedo to 0.02 < p < 0.05 (see 
e.g. Jewitt [2005]) . For 17P/Holmes an absolute magnitude, 



Table 1. Flux densities S v (250) at 250 GHz in a 11" 
beam of 17P/Holmes. A and r are the comet's distance 
from Earth and Sun at time T after the outbreak. 



date 


T 


A 


r 


Su (250) 


ref. 




days 


AU 


AU 


mjy 




Oct. 27.105 


3.3 


1.630 


2.447 


64.5±2.8 


(1) 


28.205 


4.4 


1.628 


2.451 


55.5 


2.6 


(2) 


Nov. 16.950 


24.2 


1.634 


2.530 


16.0 


4.3 


(3) 


18.769 


26.0 


1.639 


2.538 


11.6 


0.9 


(3) 


20.846 


28.0 


1.645 


2.546 


5.1 


1.1 


(3) 


23.998 


31.2 


1.657 


2.558 


5.9 


0.7 


(3) 


25.929 


33.1 


1.665 


2.566 


7.1 


1.3 


(3) 


28.142 


35.3 


1.675 


2.575 


4.5 


1.1 


(3) 


Dec. 03.838 


41.0 


1.710 


2.599 


5.5 


1.1 


(3) 


18.333 


55.5 


1.850 


2.670 


4.1 


1.5 


(3) 



(1) extrapolated from the flux of 2.3 ±0.1 mjy observed at 88.6 GHz by 
Boissier el al. (see Sect. [3l 

(2) extrapolated from 2.1 ± 0.1 mjy observed at 90.6 GHz (ibid.) 

(3) this paper 



H^, of 16.6 (Tankredi 2006}) was found and a median nu- 
clear diameter, djy, of 3.2 km derived within the limits 
of 4.6 and 2.9 km, corresponding to the albedo range. 
Meanwhile, Lamy et al. (2005) report a diameter = 
3.42 km, obtained by a single snapshot by the Hubble 
Space Telescope (HST). We prefer this direct measure- 
ment, even though it might need a correction, if the nu- 
cleus is not spherical. 

Bulk density, porosity. The bulk density may change 
from comet to comet, depending, e.g., on the outgassing 
history. For our model the value p = 0.5 g cm -3 was se- 
lected, derived by Rickman fl989) for the comet Halley 
data and from observations of 29 short period comets by 
Rickman et al. (1987). For the porosity (fraction of void 
volume/bulk volume) we assume a value 0.60. This pro- 
vides ample storage for sublimated molecular gas inside 
the nucleus. 

Equilibrium temperature. One needs to know the 
brightness temperature, T^, of the nucleus and the dust 
grains to calculate their emission. In the absence of new 
data, we assume that they both will be close to the equi- 
librium temperature, T eq . For a heliocentric distance of r 
«2.45 AU and an albedo p of 0.04, we assume T b sal 75 K, 
i.e., identical to T eq . 

3. Observations 

The outburst of comet 17P/Holmes came at an unfa- 
vorable moment, when on Pico Veleta the M Ax-Planck 
Millimeter BOlomctcr array (MAMBO) had not been in- 
stalled on the 30m telescope; in Effelsberg the 9 mm wave- 
length receiver was not operational at the 100m telescope; 
and for the Atacama Pathfinder Experiment (APEX) 
12m telescope, the comet was below the declination limit. 
About 25 days later, when MAMBO went back into oper- 
ation, we started a series of maps and ON/OFF observa- 
tions at its effective frequency of 250 GHz, trying to see the 
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aftermath of the outburst. The observing and evaluation 
procedures of MAMBO observations are standard routines 
and have been frequently reported; see e.g. Altenhoff et al. 

pro . 

The results are collected in Table 1. Prior to our mea- 
surements, the comet had already been detected with the 
Plateau de Bure Interferometer (PdBI) near 90 GHz by 
Boissier et al. (2008 , 2009}. Their results, generously made 
available to us prior to publication, are included in our 
analysis. We have scaled the 90 GHz flux densities, ob- 
tained with a synthesized beam of 5.7 x 7.3 arcsec, to the 
angular resolution of our MAMBO data (11 arcsec), and 
we extrapolated the signal to 250 GHz with the canonical 
spectral index of comets SI = 2.7, reversing the procedure 
of Jewitt and Matthews (1992) to derive the spectral index 
of comet Hale-Bopp. This method was intensively tested 
by Altenhoff et al. (|2lTO . 

Each stage of the optical development has an equiva- 
lent one at millimeter (mm) wavelengths. The optical ob- 
servations are summarized in Sekanina (|2008p . e.g. with 
the total magnitude m\ as a function of time, "the light 
curve" . 

The mm data are compiled in Table 1 and Fig. 1. The 
extrapolation of the PdBI observations to 250 GHz is fairly 
accurate, and the combined errors of extrapolation and 
observation are indicated by the size of the symbols. The 
small beam broadening by the comet, reported by Boissier 
et al. (|2008p shows that the source is optically thin. The 
two data sets are interpolated, suggesting a signal loss of 
7 % per day. The series of nuclear magnitudes TO2 shows 
a similar slope. 

In a separate paper, Altenhoff et al. (|2008p show that 
most cometary mm/radio light curves can be represented 
by the following equation: 

S v = ^, A- 2 x r" 1 ' 7 

with A and r the geocentric and heliocentric distances in 
AU, respectively. The constant S u ^= 74.5 is derived from 
the last data points. 

Thus the light curve is calculated and plotted in Fig. 
1. It is obviously a reasonable fit for the time after day 
33, when insolation and dust production (determining the 
intensity of the mm radiation) are apparently coming to 
equilibrium. For the first 30 days, this radio light curve 
is the baseline for the burst. As a further indicator of 
cometary activity, we use the nuclear magnitude, m2, re- 
ported with the astrometric positions (Marsden I2007|) . 
These values with limited accuracy are averaged over three 
days (typically over 100 observations) to reduce the noise. 
These data also confirm increased nuclear activity in the 
first 30 days. Red circles show the H2O production rates 
measured with the Solar Wind ANisotropy (SWAN) ex- 
periment on the SOlar Heliosphcric Observatory (SOHO) 
reported by Combi (2007|). This system has a beam of 
about one degree, probing the water production of about 
4 days. This may be a crude guess, considering that we 
are using observing results obtained with very different 
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Fig. 1. Comet 17P/Holmes: Compilation of spectroscopic 
and continuum observations. The black dots and the black 
diamonds represent the mm continuum data at 250 GHz, 
taken with the 30m and the PdBI, respectively. Dashed 
line: model of mm halo (see text). Red open circles: H2O 
emission, observed with SWAN, light blue dots: HCN 
emission, blue squares: HCN emission. Sec text for refer- 
ences. The dotted magenta curve shows the optical nuclear 
magnitudes mi , as an indicator of the nuclear activity. The 
spectroscopic data sets are normalized to their respective 
maximum. 

resolutions. Even though we guess that, with the result- 
ing smearing, the production rate might fit even better to 
our observed extended cometary activity. 

Spectroscopic observations of HCN by Biver et al. 
(2008) at Pico Veleta and at the Caltech Submillimeter 
Observatory (CSO), and by Drahus et al. (|2007|2008p 
with the Arizona Radio Observatory (ARO), are shown 
for comparison. The data sets are consistent which each 
other and show a steeper decay than the cometary activity 
described before. 

4. Mass determination 

Fine dust. Optically, the scattered light by small dust 
particles is dominating the appearance of comets, even 
though the mass of these particles is low. Sekanina |1982) 
has estimated the mass of 2 /xm sized fine dust in comet 
17P/Holmes (see Table 2) near its outbreak. The size of 
the scattering particles is too small to detect with radio or 
mm telescopes. This dust is responsible for the optical ap- 
pearance seen at magnitude mi . The particulate dust and 
the bulk of the molecular gas are almost invisible optically. 

Particulate dust. Radio and mm continuum observa- 
tions measure the thermal emission of dust particles of size 
> 10% of the observing wavelength, here > 0.2 mm. Since 
the observed signal is proportional to the integrated par- 
ticle cross sections, but the particle mass is proportional 
to its volume, the mass of big particles is underestimated, 
so observations at different wavelengths are needed for a 
more precise mass estimate.. We estimate the dust mass 
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with the photometric diameter to be the size of a disk 
at the distance of the comet with its equilibrium tem- 
perature, radiating as black body, yielding the same flux 
density as the radio/mm halo. For cometary dust, we find 
that the black body condition (emissivity «1) is fulfilled 
with a density of 1 g cm -3 and a layer depth of 3 wave- 
lengths (as confirmed by the rigorous halo evaluation for 
comets Hyakutake and Hale-Bopp (Altenhoff et al.[2000). 
This allows calculation of the dust mass in the halo for 
any observed signal. 

Dust production rate. The average particle moves 
through the telescope's diffraction beam in about 60 
hours, and the resulting dust production rate and the dust 
mass in the beam are listed in Table 2. 

Hypothetical pre-burst dust. The radio light curve, as 
defined above, can be extrapolated backwards over the 
whole apparition to calculate "hypothetical" signals and 
masses that would have been emitted in the absence of 
the dust cover. This total hypothetical mass is a factor 
3-4 higher than the total mass released within 33 days 
after the outburst, i.e. the burst dust mass. Thus, we can 
safely assume that the insolation provided enough energy 
to start sublimation within the nucleus. 

Accuracy estimate. Within our observing interval, the 
mass in the halo is approximately proportional to the ob- 
served flux density in the beam, thus to the observing 
accuracy. Therefore the relative accuracy from day to day 
and of the dust production rate is quite good. The abso- 
lute accuracy depends on our knowledge of the absorption 
coefficient k of cometary dust, whose uncertainty was es- 
timated by Altenhoff et al. (|2000|) to be about a factor 2. 
The accuracies of the mass determinations of the small- 
grained dust (Sekanina I2007|) and of water (Combi[2007) 
have unfortunately not been reported. 

5. Interpretation 

Start of the outburst. The nuclear structure of comets 
9P/Tempel 1 and 17P/Holmes before the outbursts are 
probably alike, a densely packed dust cover («1 m) be- 
low a layer of pure water ice (wlO m), below amorphous 
dirty H2O ice, whose upper part is possibly free of highly 
volatile molecules. At 9P/Tempel 1 the impactor acted 
as the exothermic energy source to blow off the pancake- 
shaped dust cover, as the scheme of Sekanina (|2007|) sug- 
gests, making it a type (c) split nucleus. The develop- 
ment for 17P /Holmes is different. When H2O sublimation 
started inside the porous nucleus, water vapor spread all 
over the nucleus, initiating more sublimation; deeper in- 
side, and even other molecular ices with lower sublimation 
points were heated, stored there at lower temperatures. 
The effective sublimating surface inside the nucleus, es- 
timated as excess over the emission after the burst on 
day 35 when it was near equilibrium with insolation, was 
more than 14 times the nuclear surface, corresponding 
roughly to the nuclear size of comet Hale-Bopp. The sum 
of the saturated partial pressures of all ice species was ob- 
viously breaking up the air tight dust mantle, allowing the 



cometary wind to start through the dust mantle and lifting 
dust particles, piece by piece, into the halo. The break up 
of the dust mantle is hardly spectacular, compared with 
the full start of cometary wind. 

Time scale of the outburst. Different versions exist of 
the development of the outburst. Sekanina (|2008l 120070 
refers to an explosion and a single exothermic event, 
and Biver et al. (2008]) and others report a water pro- 
duction rate, which almost ends after 3 days. In Fig. 1 
the continuum observations are plotted, showing that the 
outburst-related increased continuum emission lasted for 
about 30 days, as did the increased nuclear magnitude 
to 2 . Additional proof are the numerous photographs taken 
within the first month of the outburst; see e.g. Sekanina 
(2008J) , in which the comet appears as a filled Plerion 
rather than a shell, implying that the dust injection into 
the coma continued after the "explosion" for quite some 
time. 

Molecular production rates. Production of gas-phase 
molecules is responsible for all the cometary activity. It is 
predominantly the cometary wind of the H2O molecules, 
which lifts the dust particles from the nucleus, so a cor- 
relation between H 2 and dust production is expected. 
Usually the production of different molecules shows a fixed 
ratio, so that one can, e.g., predict the H2O production 
rate from HCN observations. Not so for comet Holmes! 
Figre 1 shows that H2O production, observed with the 
SWAN satellite, lasts at least for a month as does the en- 
hanced mm continuum emission, while e.g. the spectral 
lines of HCN, CO, NH 3 (Drahus et al. [2007] Biver et al. 
2008[ Menten, I2007j) had a big signal at the start, which 
apparently petered out dramatically after 3 days, as shown 
in Fig. 1. The reason may be the temperature/depth struc- 
ture of the nucleus, because the near surface ice might be 
free of volatile molecules. 

Mass comparison. All derived masses are collected in 
Table 2, where the total mass and the mass of the dust 
cover have been calculated with the model values. The 
mass of the 2 /jm sized dust, determined immediately af- 
ter the outburst by Sekanina (12008] ), is surprisingly high, 
compared to the total mass of the dust layer! Even the 
total particulate dust mass (grains of size 0.2 to 7 mm) is 
smaller. Integrated over the 33 days of increased cometary 
activity, the dust mass produced by the outburst is ~ 2% 
of the comet's total mass. Surprisingly low is also the H2O 
mass released in the first part of the outburst, when we 
would have expected a mass comparable to the particu- 
late dust mass. The hypothetical dust mass, calculated 
from the radio light curve backwards, is about a factor of 
3-4 higher than the total mass in the burst. The energy 
released in the burst can be provided by the insolation 
before the burst, even allowing energy losses through re- 
radiation by the dust cover. 

The total accounted mass loss during this apparition 
(mass of 2 /jm sized dust, burst dust mass, burst H2O 
mass) is in total < 3.5% of the total nuclear mass. If a 
bulk density of p = 0.5 is assumed for the outer nucleus, 
this loss corresponds to a layer of 20 m thickness which is 
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Table 2. Mass budget 



Contribution 


Mass 




Comment 




nuclear mass 


1.1 10 16 | 


J 


Model, p = 


0.5 


initial dust cover 


3.7 10 13 ! 


y 


Model, p = 


1.0 


dust (~ 2/xm) 


no 14 ! 


y 


(1) 




dust mass in halo 


4.3 10 13 ! 


y 


on day 3 




dust production rate 


2.0 10 8 J 


i/ s 


for day 3 




burst dust mass 


2.1 10 14 ! 


y 


day 3-33 




dust mass in halo 


3.0 10 12 ! 


y 


day 35 (equil.) 


H2O production rate 


3.6 10 7 ! 


y/s 


(2) 




burst H2O mass 


6.5 10 13 ! 


y 


day 1-23 





(1) Sekanina 1 120071 

(2) converted from 1.2 10 30 mols/s (Combi et al. 120071 



the same order of magnitude as found from observations 
of other comets. 

5.1. Alternative models 

Sekanina (|2007ll982p explained the outburst of comet 
17P/Holmes as a splitting nucleus, whereby the secondary 
nucleus is a fragment of a jettisoned insulation mantle of 
debris. The splitting starts with an exothermic event. His 
model considers neither particulate dust with particles > 
0.2 mm and nor the H2O production, both of which con- 
tribute at least as much mass, each seperately as doeshis 
"secondary nucleus" . 

6. Conclusion 

The historic outbursts, as discussed by Whipple (1986), 
show several similarities to the present one, suggesting 
that they happened the same way, but in 2 steps. After 
all, comet 17 P/Holmes is a comet like many others whose 
appearance is determined by sublimation of cometary ices. 
What makes it peculiar is that it had a big dust cover and 
that it seldom comes close enough to the Sun to afford a 
great display of activity. Dust covers of cometary nuclei 
are standard (see model of Horanyi et al. 1984) and do 
not indicate a splitting comet. We think that the delayed 
sublimation, as explained above, is a viable alternative 
to the theory of splitting or sudden fragmentation of the 
cometary nucleus. sectionAcknowledgement We are grate- 
ful to Dr. J. Boissier (IRAM) for communicating the 90 
GHz results to us prior to publication. We thank the di- 
rector of IRAM, Dr. P. Cox, for granting special observing 
time and the staff on Pico Veleta, Spain, for their support 
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